Phonons in random alloys: the itinerant coherent-potential approximation 
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We present the itinerant coherent-potential approximation(ICPA), an analytic, translationally 
invariant and tractable form of augmented-space-based, multiple-scattering theory in a single- 
site approximation for harmonic phonons in realistic random binary alloys with mass and force- 
constant disorder. We provide expressions for quantities needed for comparison with experimental 
structure factors such as partial and average spectral functions and derive the sum rules associated 
with them. Numerical results are presented for Ni55Pd45 and NisoPtso alloys which serve as test 
cases, the former for weak force-constant disorder and the latter for strong. We present results on 
dispersion curves and disorder-induced widths. Direct comparisons with the single-site coherent 
potential approximation(CPA) and experiment are made which provide insight into the physics 
of force-constant changes in random alloys. The CPA accounts well for the weak force-constant 
disorder case but fails for strong force-constant disorder where the ICPA succeeds. 

PACS numbers; PACS; 71.20, 71.20c 
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I. INTRODUCTION 

Much research has been carried out over the past few 
decades on the nature of elementary excitations in dis- 
ordered alloys. Many aspects of the lattice-vibrational, 
magnetic and electronic excitations in such systems have 
been intensively studied both theoretically and experi- 
mentally. Of them, the electronic problem has been cov- 
ered in most detail in recent times vifith the emergence of 
first-principles techniques which have made it possible for 
the theories to attain a much higher degree of accuracy 
and reliability. Surprisingly, this is not true for phonons 
despite their being not only conceptually the simplest 
type of elementary excitation but also the most readily 
accessible to detailed experiment. From the early 60's 
till the early 80 's there were many experimental investi- 
gations of phonons in random binary alloys |[ |[ |^, ||] 
by neutron scattering techniques. More recent experi- 
ments have been lacking, probably due to the absence of 
a reliable theory. The feature which makes the theory 
of phonon excitations difficult is the inseparability of di- 
agonal and off-diagonal disorder. The reason for this is 
that the force-constant sum rule, i.e. the force constants 
between a site i and its neighbors j obey the relation 
^ii = ~ ^ij J must be rigorously satisfied even if 
the system is disordered. In other words, a single de- 
fect at one site in the system perturbs even the diagonal 
hamiltonian on its neighbors as well, thereby imposing 
environmental disorder on the force-constants. Hence, 
any theory must include diagonal, off-diagonal, and en- 
vironmental disorder as well in order to produce reliable 
results for phonon excitations in random alloys. 

From the late 60's there were many attempts to provide 
an adequate theory of phonons in random alloys. The 
first successful, self-consistent approximation was the co- 
herent potential approximation(CPA) The CPA is 
a single-site, mean-field approximation generally capable 
of dealing only with diagonal disorder (mass disorder in 
the context of phonons). In the early 70's there were 



several studies using the CPA || which failed to es- 
tablish it as a complete answer to the phonon problem in 
random alloys. The discrepancies with experiment con- 
firmed this need for a theory which could include force- 
constant changes in addition to mass disorder. Several 
extensions of the CPA to include off-diagonal and envi- 
ronmental disorder were proposed over the next several 
years ||, |lO[ |ll|, |l^, |3[ but only in certain very special 
cases, such as the separable |^ or the additive ||l^, |ll| lim- 
its of off-diagonal and environmental disorder, were there 
successes. The more general approximations ||l2|, [l^ 
produced Green's functions which either failed to retain 
the necessary analytic properties, the translational in- 
variance of the averaged system, or were not fully self- 
consistent. Moreover, all of these extensions failed to cap- 
ture the effects of multisite or cluster scatterings which 
give rise to additional structures in quantities such as 
the spectral density functions. Later attempts which met 
with some success for real alloy systems included the re- 
cursion method |15 which can handle large clusters and 
treats all kinds of disorder on an equal footing. However, 
the recursion method is neither self-consistent nor trans- 
lationally invariant when used alone. Yussouf and Mook- 
erjce were able to provide a self-consistent general- 
ization of the CPA to include 2-site scattering using a re- 
cursion method in conjunction with the augmented space 
formalism{ASF)^^. The ASF has the proper transla- 
tional invariance, yields analytic Green's functions, and 
can handle diagonal and off-diagonal disorder. 

An alternative approach was provided by Kaplan, 
Leath, Gray and Diehl [^sj (KLGD) which is also based on 
the ASF. This approach generalized the travelling clus- 
ter approximation of Mills and Ratanavararaksha [|l9| for 
diagonal disorder to include the other kinds of disorder 
and multisite effects. Using the diagram symmetry rule 
of Mills and Ratnavararaksha and the translational sym- 
metry of the augmented-space operators, they presented 
a self-consistent multiple-scattering theory which allows 
one to work with a small number of atoms instead of 
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treating large clusters as is done in recursion. It pro- 
vides analytic, translationally-invariant approximations 
at all concentrations for diagonal, off-diagonal, and envi- 
ronmental disorder. It can be applied even to problems 
of charge transfer, lattice relaxation, and short-range or- 
der in the context of electronic excitations. However, 
they illustrated their method only with one-dimensional 
models and presented it in a very general and complex 
mathematical language. 

In this paper, we present a simple, straightforward 
formulation of the KLGD method for single-site scatter- 
ing of phonons in three-dimensional lattices and provide 
the first application of it to phonons in random alloys 
. We term it the itinerant coherent-potential approxi- 
mation or ICPA ; it maintains translational invariance, 
unitarity, and analyticity of physical properties while in- 
cluding off-diagonal and environmental disorder. In ad- 
dition to demonstrating its superiority over the single- 
site CPA and its previous extensions, we provide insight 
into the physics of force-constant disorder. Our results 
reveal the complex interplay of forces between various 
atomic species in a random environment, an important 
phenomenon which has never been addressed properly. 

In Section II we describe the theory, introducing the 
augmented-space representation and its use in construct- 
ing the self-consistent scattering theory and the single- 
site itinerant coherent-potential approximation. In Sec- 
tion III we derive expressions for important physical 
quantities such as densities of states, spectral functions, 
inelastic scattering cross-sections and their sum rules in 
terms of the configuration-averaged Green's function of 
the system. In Sections IV and V we present our re- 
sults on Ni55Pd45 and NisgPtso alloys as test cases and 
compare them with experiment. Concluding remarks are 
presented in Section VI. 



II. FORMALISM 

In this section, we briefly sketch the rationale be- 
hind augmented space, introduce its representations, and 
define the notation to be used throughout the paper. 
We present our discussions here only in the context of 
phonons. The formulation of the ICPA for other kinds of 
excitations is closely analogous. 

A. Augmented space and its representations 

The description of disordered systems conventionally 
proceeds as follows: the dynamical behavior of a system 
is described by a Hamiltonian, whereas the statistical 
behavior of the disorder is imposed from outside. The 
Hamiltonian itself does not describe the full behavior of 
the random system, but has to be augmented with the 
distribution of the set of random potentials which are 
associated with the various configurations of the system. 
The physical properties are then obtained by ensemble 



averages over configurations. The CPA and its extensions 
employ this procedure. 

An alternative procedure is that instead of looking at 
the excitations of the system as moving in a random ar- 
ray of disordered potentials, the excitations are consid- 
ered to be moving in periodic potentials in the presence 
of a 'field' which specifies the disorder. The Hamiltonian, 
expanded to include the disorder field, then by itself com- 
pletely describes the disordered system. Since the infor- 
mation on random configurations is already incorporated 
into the Hamiltonian, the configuration averaging is not 
a further process as in the mean- field approaches, but 
simply an evaluation of matrix elements. The idea of in- 
troducing a 'disorder field' to describe the random fluc- 
tuations in the system by extending the Hilbert space 
to include the disorder field and by representing the 
Hamiltonian in this new space constitutes the core of 
the augmented-space formalism. The extended Hilbert 
space which captures the random fluctuations is called 
the 'augmented space'. 

Here, we work only with a binary alloy Ac^Bc^- We 
assume that each lattice site is randomly occupied by 
an A atom or by a B atom. We wish to calculate 
the configuration-averaged values of the experimentally 
measurable physical quantities, for which we need a 
configuration-averaged Green's function. In particular, 
we shall concentrate here on the configuration-averaged 
displacement-displacement (one-phonon) Green's func- 
tion 1^ 

« it) 1 « < it) ; ut (0) », (1) 
or, after Fourier transformation to the frequency domain 

In Eqs.(l) and (2) <C ^ stands for both configuration 
and thermodynamic averaging. In Eq.(l), m,n specify 
lattice sites and a/3 the cartesian directions. u"{t) is the 
displacement operator of an atom at the lattice site n in 
the direction a at the time t. In Eq.(2) a bold symbol 
represents a matrix for which all indices are to be un- 
derstood. The semi-colon ; denotes Bose time ordering, 
m is the mass operator, $ is the force-constant operator, 
and u) is the frequency which contains a vanishingly small 
negative imaginary part. The masses are random, 

m^f = niiSafiStj, (3) 

with rrii randomly taking on the value mF if species 
r=A,B is on site i. The force-constants take on the val- 

/ aB\^^ .... 

ues [ (j)"^ j if species F is on site i and species A is on 
site j. 

It is ^ G ^ which carries all the dynamical informa- 
tions of interest, and the essential difficulty of the theory 
of phonons in random systems arises from taking the con- 
figuration average of the inverse of the matrix m.Lu'^ — $. 
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The augmented-space technique |2l|] greatly facih- 
tates this averaging. The displacements u, masses m, 
force-constants and Greens function G are defined in 
the dynamical Hilbert space ^ in which the Hamilto- 
nian of the system operates. For a binary alloy, 'J is 
augmented by the space O of all possible atomic config- 
urations of the system. The resulting augmented space 
is 



representation the occupation operators rf/^ and fjl are 
transformed to 



CA \fCACB 

\/cacb cb 
cb 



-\/cacb 

\fcACB CA 



(11) 



In O or operators are represented by symbols with 
superposed carets. In the configuration representation 
within 8, the state of site i is specified by the single-site 
state \Ai) if A is on I and by \Bi) if B is on i. With respect 
to these states, the occupation operators rj'j^ , r=A,B, 



are represented by the matrices 



IB 



= 0, 
= 



1 





1 



(4) 



(5) 



The configuration of the entire system is specified by the 
direct product of all single-site states Hi Wi)^ r=A,B. 
The mass operator for site i is given by. 



(6) 



Similarly the force-constants for sites i and j are given 
by 



In transforming from the configuration representation to 
the fluctuation representation, m' goes to m and $' to $, 
as given by Eqs.(6) and (7), respectively, with the Tf^ of 
Eq.(5) replaced by the rf of Eq.(ll). Thus the dynamical 
operators m and $ are not diagonal with respect to the 
number of fluctuations or defects in the fluctuation rep- 
resentation and can create them, destroy them, or, in the 
case of (pij, cause them to travel, or 'itinerate'. We refer 
to the movement of defects induced by the off-diagonal 
elements of the rjf as the itineration of fluctuations to 
distinguish it from the propagation of phonons. How- 
ever, these operators are translationally invariant; the 
randomness in configuration is thus captured by transla- 
tionally invariant operators in the configuration space 0. 
The rf operators constitute the disorder field referred to 
above. 

Any operator A in this augmented space can be rep- 
resented in block form. 



A A' 
A't A 



(12) 



with the cartesian indices understood. 

Consider now a rotated representation for site i in 
which the basis vectors for its configuration space are 
given by 



|0,) = V^IIA,) + 

= ^\A,)~^\B,). 



(8) 



Constructing the configuration average of any operator A 
in 8 can be carried out simply by taking the expectation 
value of A with the state 



i/>=niO'^)' 



(9) 



Thus |0i) is the site-average state (or the virtual-crystal 
state), jlj) describes a fluctuation away from the average 
state on site i, and 



i/.)=ii,>nio^-)- 



(10) 



is the state in which there is a fluctuation or a defect in 
the average state |/) only on site i. In this fluctuation 



where the bold notation A implies a matrix in the site 
and cartesian indices. The four elements of the block 
matrix are given by 



A 

A' 
A't 
A 



PAP, 

PA(l-P), 
(1 P) AP, 
(l-P)A(l-P) 



(13) 



where P, the projection operator onto the virtual-crystal 
state, is given by P = |/)(/| • Thus, we see that A is 
the configuration average of the quantity A while A',A'^ 
generate the coupling between the average and the fluc- 
tuation states and A is that part of A entirely within the 
space of fluctuation states. 

In the present paper we shall make the approx- 
imation of treating explicitly only single fluctuation 
states I/;) in the fluctuation space — |/)(/|, although 
multiple-fluctuation states are treated implicitly via a 
self-consistency condition. States in $7 can then be speci- 
fled by \ij) or \ifi) where i is the site index of the dynam- 
ical variable in 5*, position or momentum, with the carte- 
sian index understood. For the site indices of the corre- 
sponding matrix elements we shall often use the compact 
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notation 

(z/|A|j7) = 

(ifi\A\jfi') = A?'''^ 
(i/|A|j7;) = 4'\ 
(z/,|Ab7) = 4f , (14) 

where I and /' denote the locations of the concentration 
fluctuation or defect. The parentheses around I indicate 
that it is neither a site nor a cartesian direction index, 
but indicates instead the position of a fluctuation in the 
lattice. 



B. Multiple-scattering picture 

A phonon propagating in a random alloy undergoes 
irreducible multiple scattering both repeatedly off 
a single fluctuation and successively off fluctuations on 
the different sites it encounters in the process. The CPA 
takes into account the former but not the latter. To il- 
lustrate how the treatment of this process of multiple 
scattering by fluctuations differs between the CPA and 
our formalism we employ a cartoon diagram (Fig.f ). The 
top panel, a two-dimensional cross-section, illustrates the 
multiple-scattering process included in the CPA. There, 
the filled circle is a single 'fluctuation site' immersed in 
an average medium denoted by open circles. The ar- 
row on the left is the direction of phonon propagation. 
When the phonon meets the fluctuation site, it under- 
goes irreducible multiple scattering at that site. In the 
CPA(diagonal disorder), the irreducible scattering by the 
defect site is confined to the defect site. The circle around 
the fluctuation site indicates the region of influence of 
the perturbation. None of the springs are affected by 
the presence of this defect since the force-constants are 
the same everywhere. One does an averaging over all 
the possible occupations of the single site. The phonon 
diagrams of the self-energy which describe this multi- 
ple scattering process completely are shown in Fig. 2(a). 
There, the filled circles represent the fiuctuation sites, 
the dotted lines represent successive scatterings from the 
fluctuation site, and the double solid line represents the 
self-consistent propagator. The lower three panels in 
Fig.l illustrate scattering sites in the ICPA. The differ- 
ence from the CPA is that the region of influence is not 
only the site of fluctuation but also its neighboring en- 
vironment around the fluctuation site. The flgure shows 
an example (dotted contour) where the environment in- 
cludes nearest neighbors only (The calculations could be 
extended to further neighbors as well). When the phonon 
interacts with the fluctuation site in the top panel of the 
three, it scatters also from all of its neighbors since their 
spring constants also undergo changes (denoted by the 
thick spring lines in contrast to the thin ones for the av- 
erage medium). The whole cluster of atoms undergoes 
fluctuations in force-constants as the occupation of the 



Multiple Scattering Picture 



CPA(mass disorder only): 




ICPA(Mass and force constant disorder): 




FIG. 1: Multiple scattering picture in the CPA(top) and with 
the ICPA(bottom). The filled circle is the site of the fluctua- 
tion, and the contours around it indicate its area of influence. 
The arrows with the ICPA indicate the itineration of the fiuc- 
tuation to neighboring sites. The details are given in the text 
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(a) CPA 



b) ICPA 



FIG. 2: The self-energy diagrams in the CPA(top) and a typ- 
ical example in the ICPA(bottom). The details are given in 
the text . 



fluctuation site changes. One has to keep in mind that 
the force-constant between the fluctuation site and its 
neighbor on the right, say, depends on the occupation of 
both sites, as is true for the next neighboring site on its 
right as weU. So, one is led to include the irreducible scat- 
terings by the fluctuation on all neighboring sites, which 
then requires inclusion of scattering by the fluctuations 
on its neighbors etc., until the irreducible scatterings ex- 
tend throughout the entire sample. A simple example of 
this process is indicated in the middle and bottom panels. 
Indeed, Mills and Ratanavararaksha[|^ have shown that 
once there are non-diagonal terms in the scattering, the 
self-energy must include these migrations (itinerations) 
of the scatterer throughout the sample in order to attain 
unitarity and thereby guarantee that the average Green's 
function will be properly analytic or Herglotz. The self- 
consistent scattering and the resulting coherent potential 
about a single defect thus itinerates from defect to defect 
throughout the sample, making it an itinerant coherent 
potential. The scattering could have started from any 
site in the sample so that the result is also fully transla- 
tionally invariant, and the self-energy is g-dependent but 
diagonal in the g-space of the Brillouin zone of the un- 
derlying periodic lattice structure. Fig. 2(b) illustrates 
a typical self-energy diagram in the ICPA. The solid and 
dotted overlapping ellipses denote the multiple scatter- 
ing by a single-site and its neighbors, i.e. by a cluster 
of atoms, and the subsequent itineration of this process. 
The thin dotted lines and the thick double lines are as in 
Fig.2(a). 

In the multiple scattering framework, we calculate the 
self-energy E (ui^) , defined by 

« G (c.^) »= [G-i - S {u;')]-\ (15) 

where Guca is the unperturbed Green's function, 

(16) 



Gyca = (mw^ _ $) 



and m and $ are the configuration-averaged mass and 
force-constant operators respectively. 



Our major task is to calculate the self-energy S (i^^). 

Let us consider K = ^mw^ — <I>^ = G~^. Using the 

2x2 block representation of augmented-space operators 
of Eq.(12) we get 



G = 



( ^ 


G'^ 






[ G't 


G J 




K 1 



(17) 



Using the relation for the inverse of an operator in 2 x 2 
block form [E6[, namely. 



(a - A'A-iA't 
- fAA't^^A- A 



AA'^^A 



A't 



A- A'tA-iA' 



(18) 



we get 

(riic^2_5^]|-iK't]-i 
= [G-i-K'FK't]-i. (19) 
Therefore, the self-energy is given by 

S = K'FK't, (20) 

where 



K 



{G-ll-Y}-\ 



and where 



nil — m I w 



(01 - 



(21) 



(22) 



The quantity V denotes all perturbations to the average 
medium, and F contains the itineration of the fluctuation 
in the average medium. 

Up to this point, the scattering formalism is exact. 
We now introduce the ICPA by restricting the states 
within the configuration space O — |/)(/| to the single- 
fluctuation states, the notation for which is given by 
Eq.(14). Making the site and cartesian indices explicit, 
we obtain for E in Eq.(20), under this restriction, 



_ \ " T^{m) p(m)(n) 7^/t(n) 



(23) 



The sumations are over the repeated indices, and the 
fluctuation itinerator F is given by a Dyson equation. 



fW(^-) = G_[<5(,)(,)+^vW«f(')W], 



(24) 



where only the site index of the fluctuation is shown. The 
quantities in (23) are translationally invariant as follows: 



ik i — m.k—m'> 



F, 



(m)(n) 



kl 



k—mA—m ' 



(25) 
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The single fluctuation in Eq.(23) can be considered to 
have been 'created' by K'*^"^ at site n, itinerated to site 
m by F(")(") and 'destroyed' by K'^K™) at site m. The 
K, K'''' and F matrices have elements which are non-zero 
only for site indices within the environment of the ap- 
propriate defects i.e. the indices i and k (l and j) must 
be within the neighborhood perturbed by the defect at 
m{n). The terms with more than one fluctuation(defect) 
present at a time correspond to coherent pair and 'de- 
fect cluster' scattering and are neglected in the single-site 
scattering considered here. All of these operators act in 
the augmented space. The Equations (20)-(24) define an 
itinerant single-site multiple scattering theory. 



C. Self-consistency 

The restriction in Eq.(23) to states of — |/)(/| con- 
taining only a single fluctuation is a very severe ap- 
proximation. Multiple-fluctuation states are of course 
present in F and contribute to E. In the spirit of the 
CPA, these are included approximately by introducing 
self-consistency. As in the CPA ^ we obtain self- 
consistency by replacing Gvca in F in Eq. (24) by a condi- 
tional propagator G^'\ identical to G =<^ G 3> except 
that all irreducible scatterings beginning or ending on 
site i are omitted, so that F would then be given by 

fWO) = G(*)[%)0) -f ^ V««F(')(^)]. (26) 

In parallel with Eq.(15), G*^*) contains a conditional self- 
energy which, is like Eq.(23), except that it includes 
only those scatterings that neither start nor end on i, 

G« = [(G„,a)-^-S«]-\ (27) 

^(j) ^ ^ x'^'^F^'^^^^K't^™). (28) 

Referring to Fig. 2a, the double line in the multiple- 
scattering graphs is the propagator G'*-' when the solid 
dot refers to site i. We obtain the itinerant CPA by al- 
lowing K',K'''', and V to include force-constant disorder 
as well and therefore defect itineration in Equations (26)- 
(28). This closed set of equations defines our single-site, 
self-consistent, multiple-scattering theory which, when 
solved, yields F. Inserting F into Eq.(20) for E and the 
result into Eq.(15) then yields G. It is already known 
that Eqs.(15),(20),(26)-(28) have a unique solution which 
yields a Herglotz average Green's function [Q. A ma- 
jor difference between this and previous generalizations 
of the CPA is that for scattering from single-site fluc- 
tuations with off-diagonal and/or environmental disor- 
der, as is considered here, the matrix representation of 
the operator V has elements which transfer or itinerate 
the fluctuation from site to site. This feature causes the 



self-energy to have nonzero off-diagonal elements in real- 
space extending across the sample and thus contributes 
importantly to such quantities as the two-particle vertex 
corrections in a way the CPA cannot |^ . 

It now remains to solve these equations, making use 
of the translational symmetry of the augmented-space 
operators. We accomplish this by Fourier transforms on 
the fluctuation-site labels, 

A {q)rnn = E ' (29) 

1,1' 

and 

where Rw is the lattice vector connecting the fluctuation 
sites I and I' , m and n are neighbors of I and I' respec- 
tively, and the q sum is over the Brillouin zone. 

We can also effect Fourier transforms on the site indices 
themselves. That of the self-energy is 

E(g)=iV-i^E,,e-^"'-«--. (31) 
From Eqs.(20), (25) and (29),it follows that 

In this notation, Eq (26) becomes 

F (<7l™„ = G(oi + E G<^ly (9),.p F {q)^„ . (33) 

rp 

The cartesian indices here are implicit so that each quan- 
tity is a 3 X 3 matrix . 

Since the range of interaction in real-space is finite, the 
perturbations K'^'^ and V (g) are finite matrices, nonzero 
only over a finite set of real sites. For example, if we 
consider nearest-neighbor perturbation only in a single- 
site approximation, the V {q) and K' are 3(Z -I- 1) x 3 
(Z + I) matrices where Z is the number of nearest neigh- 
bors. This is the minimum matrix size necessary to ex- 
hibit all the impurity modes or states about each fluctu- 
ation site. 

In full matrix notation, we obtain 

F{q) = [G^"^-' -V{q)]-K (34) 

These matrices , for example, for an fee lattice, are of 
dimension 39x39. 

In order to evaluate G^"-* we rewrite Eq. (27) as, 

G(") = [(G.,,)-^ - E(")]-i = [« G »-i +E(°)]-i, 

(35) 

where, E^'^' = (E — E'^*''). The conditional self-energy 
E'-'') contains only those scatterings which either start or 



7 



end with a perturbation caused by a fluctuation at site 
0. Thus, to evaluate the self-consistent propagator G^"\ 
we need to know ^ G 3>. But <C G ^ is obtained from 
Eq.(15), which becomes 

«G(g)» = [G,ca{q)-^ -^m'\ 
< Gy > = iV-i ^ < G (g) > e-'«-^- . (36) 
q 

After reaching self-consistency by the procedure de- 
scribed below, we use these expressions to calculate den- 
sities of states and spectral functions. 

The conditional self-energy £(0) can be broken up into 
two contributions: 

(i) Scattering that starts from a defect at site and 
ends at site j. 

(ii) Scattering that starts at j but ends at 0. 
This decomposition results in 

_K'(o)f("^(°)K''^(°). (37) 

The last term is subtracted to avoid overcounting when 
J=0. 

In a block notation similar to that of Eq.(12), we have 



;(o) 



Si E3 
A 



G 



(0) _ r g1«) G 



i(0) 



ptio) MO) 



«G»=( ^1 



(38) 



(39) 



(40) 



where, for a genaral operator A, Ai begins and ends with 
scattering about site , A2 neither begins nor ends with 
scattering about site and A3 (A3) begins (ends) with 
scattering at the site and ends (begins) with scattering 
about a site different from 0. The term S2 is since E*^"^ 
must begin or end at the site 0. From Eq. (35), we have 



g'*') =<g> (i + e^'^) <G> 

which leads to 



(41) 



G 



where 



(0) 



X[I + El - E3G2ET X + E3GT]-i, (42) 



(43) 



X= (I+G3EJ 



Gi. 



after a lengthy algebraic analysis which was previously 
given in Ref.18. 



In order to evaluate these expressions [Q, we need to 
calculate four terms: Gi,Ei, G3E3 and E3G2E3. The 

" ' ' t a. finitp «nm nf finitp ma.trir'PS a.nd 



calculate tour terms: iori,Zji, tjrsZjg _j„^_3. 
first term Si is just a finite sum of finite matrices and 
can be evaluated directly, but the other two terms involve 
sums which range over all sites in the solid and must be 
evaluated by Fourier transforms. This is done in the 
following way, 



(G3Et)^_^^ - E E « G »t..m^r K't^ 

m r,n,l 



«i^r-E<G(.^) >,.e!3, 

r 



which becomes 



a r 



(44) 



and, similarly, 
E3G2E3 



t,t' 



q,r,r' 

-(EiGiEi + EiG3Et + E3GjEi)i,ts (45) 



where 



M 



^ 1 



(46) 



In these equations, is the index of the single fluctuation- 
site in consideration; r,r' ,t,t' ,l,n are the neighboring 
sites of 0; and m, m' are general sites in the sample. So, 
it is clear that one needs to work only on matrices of 
size 3{Z -I- 1) x3 {Z + 1) and use the Fourier transform 
of operators to handle the itineration of the fluctuation 
throughout the entire sample. An interesting point to 
note is that the quantities G3E3 and E3G2E3 represent 
the scattering and itineration of the disturbance includ- 
ing the effect of the off-diagonal and environmental disor- 
der. In case of diagonal-disorder only, they vanish giving 
G^"^ = Gi (I + Ei)"\ which is the CPA self-consistent 
propagator, and the self-consistent set of equations re- 
duces to the CPA equations. 

The inputs to the self-consistency cycle are G^"^^^ = 
G„ca(or some better guess) , K',K''I' and V(g). The 
procedures for evaluating the latter three quantities are 
given in the Appendix. The cycle consists of the following 
steps : 

1. calculation of F (q) using Eqs.(33) and (34). 

2. calculation of E (q) using Eq.(32). 

3. calculation of ^ G ((7) ^ and <C G (u^) ^ using 
Eq.(36). 

4. calculation of G^"^ using Eqs.(42),(43),(44) and 
(45) 
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5. If the results of steps 1.-4. are acceptably close to 
those of the previous cycle, stop. If not, use as 
input to step 1 and iterate. 

The iterations are done till self-consistency is achieved 
for each g-point in the Brillouin zone. In the process of 
achieving self-consistency, one calculates <C G in both 
real-space and in g-space; each is needed to obtain den- 
sities of states and spectral densities respectively. In the 
next section, we describe how these are used to calculate 
physical quantities of interest and discuss their signifi- 
cance. 

III. IMPORTANT QUANTITIES; SUM RULES 

In this section we derive rcsiilts for important physi- 
cal quantities such as the densities of states (partial and 
total), spectral densities (partial and total) and inelas- 
tic scattering cross sections (coherent and incoherent). 
All of these are derivable from the real-space and the 
g-spacc configuration-averaged Green's function and en- 
able us to make direct comparisons with experimental 
measurements. 

A. Densities of states 

The total density of states for a 3-dimensional system 
is defined as, 

^ (w) = ^^"^{^^ « mG(a;2) »}, (47) 

where m is the mass matrix, and N is the number of 
sites. In augmented space wc have, 

<mG>ii = (i/|mG|z/), 

= (z/|m|z/)(z/|GK/) 

+(z/|m|z/,)(z/,|G|z/), 
= mGoo+m'(i/i|G|i/). (48) 

To evaluate the second term, we use the notation of 
Eq.(12) for the operators G and K = G~^. Then, using 
Eq.(18), we obtain 

(i/i|G|i/) = G't = -K-^K'tG = -FK'tG. (49) 
We can, therefore, write 

I j,n 

Fourier transforming over the fluctuation site according 
to (29) gives 

mG\if) = -^T.T.Pi^o,j-ie''-^" 

l,j,n q 



The Fourier transform of <C G on the real-site index 

now gives, 

X « G(gO » e^?-«" 

Finally we obtain, 

mG\if} = -EE^®o,,n^V^'"'^''«G(9)»' 

mp q 

where, m = j — I, p = n — I, the neighboring sites per- 
turbed by the fluctuation. All the terms on the right 
hand side have been calculated already in the process of 
achieving self-consistency. The evaluation of the average 
density of states is thus straightforward. 

The partial density of states for atoms of type s is given 

by 

771 

y (^). = ^^^{T^^ « G{^''r »»}, (50) 

where 

« G'' G'' >oo= (0/|^^G|0/) (51) 

because of translation invariance. Wc thus have 

« »o=« G^^ »oo = (0/|^^|0/)(0/|G|0/) 

+ (*/|57g|0/o)(0/o|G|0/)(52) 

and, from Eq. (11), it follows that 

^Ma = -^Jm[cA{(0/|G|0/)} 

+y^{(0/o|G|0/)}], 
^Ms = -^Xm[cB{(0/|G|0/)} 

-x/^I^{(0/o|G|0/)}]. 

(53) 

The elements of G in Eq.(53) were already evaluated 
while calculating the average density of states above. 

The partial Green's functions <S G" »o are used in cal- 
culating the incoherent scattering structure factor which 
is directly measured in the experiments, 

« Sin^oh (Q, »= E -Im^G' (w) >o -Q, 

(54) 

where bs is the incoherent scattering length for atoms of 
type s, and Q is the phonon wavenumber. 

B. Spectral densities 

The average spectral function is defined as, 

< Ax {q, w^) >= -Im «: Ga {q, w^) >, (55) 
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where A is a normal-mode branch index. More interest- 
ing quantities to calculate are the conditional or partial 
Green's functions <C G*' ((f, w^) » in g-space because 
these enable one to calculate the coherent-scattering 
structure factors which are measured directly in the 
neutron-scattering experiments and are given by, 

< Sx {q,(jj) :S>coh^^dsds'-Im < Cf {q,oj^) >, 

ss' 

(56) 

where, ds is the coherent scattering length for the species 
s. 

The conditional Green's functions are defined as 
« G''' {q,uj^) »= ;^ ^ « G''' (to^) e-«-^«. 



= {^f\vm)mG\Jf){Jf\fif\Jf) 

+{^m!mmG\jf){jf\rf;\jf) 
+{mvn^.f){mG\jf,){jmf\jf) 
+{^mtm{mG\jf,){jf,\rjf\jn57) 

In Eq.(57) the index A is to be understood. These four 
terms include all the possible scattering processes when 
two different sites arc occupied by two species. The four 
different terms involve calculations of the Green's func- 
tion under various circumstances of coupling between the 
average and the fluctuation states weighted by the appro- 
priate concentrations . 
We obtain from (57) 



where 



J2 Fon {q) K'^me''-^- « G {q,u;') », 

nm 

^ « G (q, iv') » e-^«-^"<„F„o (g) , 

nm 

Foo (9) + E E ^on {q) K'^me''-^- 

nm Ip 

X « G {q,w^) » e-'^^-^'Kl^Fpo {q) , (60) 



and where n, m, I and p are the neighboring sites of the 
fluctuation site influenced by the perturbation. For 
the lattices with each site having inversion symmetry, 
Ti=T2 holds because Ti and T2 are the contributions 
from two processes which are conjugate to one another. 
In that case, «; G^^ (9, w^) >= « G^^ (9, w^) » when 
ca=cb- 

The sum rules for the conditional Green's functions are 

derived the following way: Integrate oj G {^^), Eq.(2), 
along the real axis, closing the contour above at infinity, 
obtaining 



dii) u G (w^) = m ^ni, 



or 



CsCs'{if\G\jf) + [Cs'\/Cs{l - Cs) 



/ doj 2<jj ImG (oj^) = m-^TT. (61) 
Jo 

Similarly, using Eq.(57), taking the Fourier transform of 
Eq.(58), carrying out the contour integral and inserting 
(61) yields the sum rule for the partial spectral functions 

f°° / c 

/ du> 2uj Im <t:G'^ {q,(j^) :$>=n—Sss'. (62) 



X (-1) 



{ifi\G\jf)] + [cs\/cs'{l — Cs') For the total Green's function we obtain 



X (-l)(l-"°)(^/|G|J7,)]+' 



x(-l) 



:\h\G\j.f, 



(58) 



The integer is equal to 1 if s=A and is equal to if 
s=B. 

These terms can be easily calculated using Fourier 
transforms as has been previously demonstrated for the 
density of states. The final forms of the conditional 
Green's functions in g-space are 

< G^^ {q,uj'^) > = 4 < G (g, u"^) » +CAy/^X^{Ti + 

T2) + cacbT^, 

< G^^ {q,LO^) » = 4 < G {q,LO^) > -cbVca^{Ti + 

T2) + CaCbTs, 

«; G^^ {q,uj'^) > = CaCb «: G {q,oj'^) > +y/EX^{cBTi 

-CAT2) - cacbTs, 

< G^^ {q,uj'^) > = CACB < G (g,u;2) > -^/c^{caTi 

-CBT2) - cacbTs, (59) 



r 2u Im < G (q,u^) >= tt (— + —] . 
Jo KmA mB J 

(63) 

The experimental dispersion curves are obtained from 
the wave-vector dependence of the peak frequencies of 
the structure factors as measured, after a deconvolution 
of the experimental resolution function. The question is 
whether the dispersion curves so obtained, which incor- 
porate the effect of the coherent scattering lengths, differ 
significantly from those obtained from the peak frequen- 
cies of the Green's function itself which gives, in princi- 
ple, a proper description of the dynamics but docs not 
contain the scattering length weighting. To answer that 
question one needs to recognize that the peak positions 
in ImG are very closely related to the zeroes of TZeG^^ 
at a given wave vector. If we diagonalize the Hermitian 
TZeG~^ both with respect to mode and species index, 
each of the two species components of TZeG~^ will have 
a zero. Correspondingly, each of the two components of 



ImGs will have a peak, if XmS does not wipe it out. 
So, in the species representation, the different matrix el- 
ements XmG""* will thus all have these peaks at nearly 
the same frequencies. Thus, the weighting of the G^^ 
by the scattering lengths will not shift the peak positions 
significantly even when the scattering lengths differ ap- 
preciably. The intensities and the lineshapes of Scoh and 
XmG may differ significantly, but the peak positions will 
generally differ little. In summary, the structure factor 
fairly accurately reflects the phonon dynamics contained 
in G with regard to the dispersion curves, an important 
fact illustrated below in the next two sections where we 
present our calculations on Ni55Pd45 and Ni5oPt5o alloys. 



IV. APPLICATION TO NI55PD45; WEAK 
FORCE-CONSTANT DISORDER 

In the next two sections we explore the relative impor- 
tance of force-constant disorder and mass disorder for the 
vibrational properties of random alloys in two specific al- 
loys Ni55Pd45 and NisoPtsg. In the former alloy, the mass 
disorder is much larger than the force-constant disorder. 
The mass ratio mpd/rriNi is 1.812, whereas the Pd force- 
constants are only about 15% larger than those of Ni 
p2| , |2^ . In the latter alloy, both the mass disorder and 
the force-constant disorder are large, providing an inter- 
esting contrast between the two materials. For both the 
cases we have done our calculations on 200 w-points and 
have used a small imaginary frequency part of -0.01 in 
the Green's function. For the Brillouin zone integration 
356 (f-points in the irreducible 1 /48-th of the zone pro- 
duced well converged results. The simplest linear-mixing 
scheme was used to accelerate the convergence. For both 
the cases the number of iterations ranged from 3 to 13 
depending on the frequency w. 

For Ni55Pd45, we compare the results of virtual crys- 
tal(VCA), CPA, and ICPA computations, using the 
ICPA force-constants to construct the averages used in 
the VGA and GPA and compare the results with exper- 
iment. We make a distinction between that use of the 
VGA and of 'mean crystal' models in which the aver- 
age mass is employed and a set of 'mean-crystal' force- 
constants are fitted as parameters to the experimental 
data. 

Kamitakahara and Brockhouse ||^ investigated 
Ni55Pd45 by inelastic neutron scattering and reported 
a strange observation. A theoretical calulation based 
on a mean crystal model having the average mass and 
fitted force-constants between those of Ni and Pd agreed 
closely with the experimental dispersion curves. This 
was quite a puzzle because it suggested that the large 
mass disorder had little effect. There were theoretical 
studies on this system using recursion ||2^ and the 
average t-matrix approximation ]25[| , but no theoretical 
results for the frequencies were available. In an attempt 
to solve this puzzle, we have carried out calculations 
with the GPA, the IGPA, and the VGA, as well as with 



10 

[^00] [^^0] [t^t^Q 



9 



6 
3 



6 







L 06' 



</c/ T 

/, , , 


— ^ \n X) - Q 

w 


y 

9-/ - 

p'/ 
,/ 

/ yJS ° - 


. . . 


T, 


pe< 
/ y 

p/ 
/ 

1 ^<> ° - 



1 ?- 0.5 



(b) [^00] [m vm 




FIG. 3: (a) (Top panel) Dispersion curves (frequency u 
vs. reduced wave-vector Q for Ni55Pd45 calculated in the 
ICPA(soIid line) and in the VCA(dashed line). The circles 
are the experimental data [^. (Bottom panel) Dispersion 
curves for Ni55Pd45 calculated in the CPA(solid line) and in 
the VCA(dashed line). The force-constants used are given in 
the text. The circles are the experimental data |3|. (b) Dis- 
persion curves for Ni55Pd45 calculated in the CPA(solid line) 
and in the mean-crystal model(dashed line) using the force- 
constants of Ref. 3. The circles are the experimental data 

i- 



the mean-crystal model used in Ref. 3 . 

For the IGPA calculation, we assumed that the ex- 
plicit scattering caused by the force-constant disorder 
was confined to the nearest-neighbors. This assumption 
is justified because the nearest-neighbor force-constants 
are an order of magnitude larger than those of the fur- 
ther neighbors so that the nearest-neighbors feel the ef- 
fect of disorder most strongly. For the virtual-crystal or 
the average medium into which the scattering was em- 
bedded, we kept terms in the Hamiltonian up to the 
fourth neighbor, which turned out to be sufficient. The 
problem with force-constant disorder scattering calcula- 
tions is the general absence of prior information about 
species-dependent force-constants. We note that Pd is 
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the larger atom here. In the alloy, the Ni-Pd separa- 
tion is larger than the Ni-Ni separation. As a result, 
the Ni-Pd force-constants should be less than the Ni-Ni 
ones. Using this intuitive argument and for simplicity 
in this illustration, we kept the </'^'^_^^ and 4>'pd-Pd 
same as those of the pure materials |^ and reduced the 
't''m-Pd below the 4>N^i-Ni an a/? independent factor. 
The dispersion curves were obtained from the ICPA cal- 
culations using (j)^._p^=0.7(j)^^_fq^ (solid lines) for the 
nearest neighbors and using the force-constants of Ref.3 
for the higher neighbors. They are compared in the top 
panel of Fig. 3(a) with the experimental results ||] (open 
circles) and the VGA for the same force-constants(dashed 
lines). These ICPA dispersion curves were constructed 
by numerically determining the peaks in the coherent 
scattering structure factor <C Sx {q, uj) ^coh given by 
Eq.(56), which was calculated using the partial spectral 
functions of Eqs.(59) and (60) and weighting them with 
the coherent scattering lengths for Ni and Pd. We could 
thus make a direct comparison with the experimental re- 
sults because the neutron data observed in the experi- 
ments inherently incorporates the effect of the scatter- 
ing lengths of the species. Excellent agreement of the 
ICPA with experiment was obtained for all three sym- 
metry directions and for each branch by varying only 
one parameter in the force-constant matrix. This sug- 
gests that the force-constant disorder is weak and the 
system is dominated by the mass-disorder, as one would 
expect from the numerical values of the parameters. It is 
confirmed by the results of the CPA calculations shown 
in the bottom panel of the Fig. 3(a), using the same 
force-constants. As in the top panel, the solid lines 
are the CPA results, the circles are the experimental 
points, and the dashed lines are the VCA results. The 
agreement with the experiment again suggests the dom- 
inance of the mass-disorder, but there are more inter- 
esting points to note. In the long-wavelength (low q) 
regime, the VCA, the CPA, and the ICPA curves are in- 
distinguishable because the self-averaging of both mass 
and force-constants over a single wavelength reduces both 
the CPA and the ICPA to the VCA. But, as we move to 
high wave-vectors , the VCA deviates to frequencies be- 
low the experimentally observed ones. This fact is due 
to the use of an average mass in the Hamiltonian. In 
the high-wave-vector region, the lighter atoms, i.e. Ni in 
this case, dominate and push the frequencies up. That 
is why the CPA and the ICPA agree very well across the 
Brillouin zone while the VCA fails for the high wave- 
vectors. The reason that Kamitakahara and Brockhouse 
got a very good fit to the experimental points in Ref. 
3 by using their mean-crystal model is that they ob- 
tained parametrized force-constants which were higher 
than those calculated in the VCA. Though they had used 
the average mass in their calculations the higher values 
of the force-constants (They used $ = ca^aa + cb^bb 
rather than $ = c\j^^aa + c%g^BB + 2cacb^ab) com- 
pensated for their omission of the effect of the mass fluc- 
tuations. This is illustrated in Fig.3(b). There, the 
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FIG. 4: (a)Partial and total spectral functions calculated 
in the ICPA for various C, values in the [C,0,0] direction in 
Ni55Pd45. (b) Partial and total structure factors calculated 
in the ICPA for various C, values in the [C,0,0] direction in 
Ni55Pd45. The solid lines are the total contribution, the dot- 
ted lines are the Ni-Ni spectra, the long-dashed lines are the 
Pd-Pd spectra and the dot-dashed lines are the Ni-Pd con- 
tributions. The details are given in the text. The left(right) 
column is for longitudinal(transverse) modes. 

dashed lines are their mean-crystal model calculations, 
the circles are the experimental points, and the solid lines 
represent a CPA calculation with the force-constants 
used in Ref. 3. Here we see that the CPA yields frequen- 
cies that are too high in the large wave-vector region. 
The CPA captures the effect of mass-fluctuation and the 
domination of the Ni atoms for higher wave-vectors, but 
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the higher values of the assumed mean force-constants 
pushes the frequencies up further, thereby worsening the 
agreement with the experiment. But, there was yet an- 
other puzzle to solve. Another striking feature is that 
in spite of incorporating the scattering lengths in our 
calculations in determining the frequencies there was lit- 
tle change in the ICPA results with respect to the CPA 
results even though the coherent scattering lengths of 
Ni and Pd differ significantly( The coherent scattering 
length for Ni is 1.03 while that of Pd is 0.6). This lack 
of change can be understood from a comparison between 
the partial and total spectral functions (Fig. 4(a)) and 
the partial and the total coherent structure factors (Fig 
4(b)). In these figures, we have shown examples of ICPA 
spectral functions and structure factors along the [C,0,0] 
direction, C = | , for a low C , one in the middle of 
the Brillouin zone and one at the edge of the Brilluoin 
zone. In Fig. 5(a), the solid lines are the total spec- 
tral function while the dotted lines are the Ni-Ni spectra, 
the long-dashed lines are the Pd-Pd spectra and the dot- 
dashed lines are the Ni-Pd contribution. In each case, 
the peaks corresponding to the dominating species and 
that in the Ni-Pd curves occur at the same general po- 
sitions. For example, in the [.3,0,0]-L curves, the peak 
in the spectral function is mostly that of Pd atoms while 
for the [.5,0,0]-L and [1,0,0]-L curves, the contributions 
are from Ni atoms, the Pd-Pd contribution here is much 
less and that too is almost completely neutralised by the 
Ni-Pd contribution in the low frequency region. The 
occurence of the peaks of the Ni-Pd spectral functions 
and that of the Ni-Ni or Pd-Pd spectral functions almost 
at the same position across the Brillouin zone suggests 
that the inclusion of scattering lengths would primarily 
alter the relative weights of various contributions and 
thereby altering only the line shapes, while the disper- 
sion curves wouldn't undergo significant change. This is 
demonstrated in Fig. 4(b). There, the solid lines are 
the coherent structure factors <C S\{q,Lu) the dot- 
ted lines are the Ni-Ni spectral functions weighted us- 
ing the scattering length of Ni, the long-dashed lines are 
the Pd-Pd spectral functions weighted using the scatter- 
ing length of Pd and the dot-dashed lines are the Ni-Pd 
spectral functions weighted by the scattering lengths of 
Ni and Pd according to Eq.(56). One can see that the 
weighting affects primarily the peak heights. These ex- 
plicit numerical results confirm the qualitative argument 
given at the end of section III. 

The disorder-induced widths are important because 
the effect of disorder is often manifested in them more di- 
rectly than in the frequencies. Kamitakahara and Brock- 
house extracted full widths at half maxima(FWHM) from 
their neutron groups by assuming that the observed line 
shape could be adequately approximated by the convo- 
lution of a Gaussian resolution function(representing the 
experimental resolution) with a Lorentzian natural line 
shape. Thus, for a comparison of our results with theirs, 
we have fitted our structure factors to a Lorentzian to 
extract the widths. The results are shown in Fig. 5. The 
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FIG. 5: Disorder-induced FWHM's in Ni55Pd45 calculated 
in the ICPA(solid line) and in the CPA(dotted line) using 
the force-constants of Fig. 3(a). The circles are the widths 
extracted from the experimental results Q. The filled circles 
had better experimental resolution. 



solid fines are the widths (FWHM) obtained in the ICPA, 
and the dotted lines are those obtained in the CPA. 
The circles are the experimental points, the filled circles 
having better resolution ^ were those in which Kami- 
takahara and Brockhouse had more confidence. Gen- 
erally, there is little difference between the widths ob- 
tained in the CPA and the ones obtained in the ICPA. 
In all three symmetry directions and for all branches, 
the ICPA performs slightly better than the CPA for high 
wave-vectors. The worst agreement with the experiment 
is for high wave-vectors in the [C,0, 0] and C] longi- 
tudinal branches and the [C,^, ^] transverse branch. In 
these cases, the low values of the widths in the theoreti- 
cal calculations can be understood from the shape of the 
structure factors. From the examples in Fig. 4 one can 
see that the agreement with experiment is good when 
we have a symmetric line shape, for example, for the 
[.5,0,0]L mode. On the other hand, the worst agreements 
with the experimental widths are for cases where we ob- 
tain a highly asymmetric line shape, for example, for 
the [1,0, 0]L mode. Fitting Lorcntzians to such asymmet- 
ric line shapes is not condusive to meaningful values of 
the FWHMs. Because they obtain higher widths than 
the theories in those particular cases where the observa- 
tions have been made with worse resolution(open circles 
in Fig. 5), it is also not clear how trustworthy their treat- 
ment of the resolution function is. 

The discussion above clearly tells us that for Ni55Pd45, 
the dominant effect is mass disorder. That alloy therefore 
does not provide a proper test of the ICPA. Nevertheless, 
our discussions show how a mean-crystal model can com- 
pensate for the neglect of mass fiuctuations in alloys with 
little force-constant disorder through the introduction of 
erroneous mean force-constants, a classic case of cancel- 
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lation of errors. 



V. NI50PT50; STRONG MASS AND 
FORCE-CONSTANT DISORDER 

The niass ratio mpt/m^i is 3, quite large compared 
to that in the NiPd system. The force-constants of Pt 
are, on an average, 55% larger than those of Ni. This 
makes it a potential example of strong force-constant dis- 
order. Tsunoda et a/.[Q investigated Nij-Pti-^, by inelas- 
tic neutron scattering and compared their observations 
with the CPA. Here, for illustration, we have considered 
x=0.5 only because that makes it a concentrated alloy 
and the failure of CPA was, qualitatively, very prominent 
at this concentration. They compared their incoherent 
scattering data with that of the CPA which predicted a 
split band separating out Ni and Pt contributions with 
a gap between them. The experiments did not reveal a 
split-band, and it was very clear that the inter-species 
forces play a significant role. We performed calculations 
with the CPA and the ICPA. As before, we used the 
ICPA force-constants in the CPA. The choice of ICPA 
force-constants was more difficult than for the NiPd be- 
cause of the larger size difference between Ni and Pt. In 
this alloy, the Ni-Pt separation is also larger than the 
Ni-Ni separation. As a result, the Ni-Pt force-constants 
should also be less than those of Ni-Ni. Moreover, a 
pair of Ni atoms would find themselves in a cage partly 
made of larger Pt atoms which would therefore reduce 
the Ni-Ni force-constants relative to their values in the 
pure material. Similarly, the bigger Pt atoms find them- 
selves compressed between much smaller Ni atoms, which 
would increase the Pt-Pt force-constants with respect to 
their values in pure Pt. Using this intuitive argument, we 
found that the following guesses for the force-constants 
worked weh: 4''j^^_j^^, (l^^'pt-pv 't^Ni-Ni and 0S-pt are 
kept the same as those of the pure materials p2i and 

4'^Ni-Ni = ^m~Ni — ^■^4'm{pure)^ 
't'^Pt-Pt ^ 4'n-Pt = l-l'/'f''t(p«re)' 

(I'm-pt = ^■^'Pm-Ni^ foi' all a,/3. 

In Fig. 6, we compare the ICPA results for the incoher- 
ent neutron structure factor (Eq.(54)) with those of the 
CPA and the experiment Q|. In Fig. 6(a), the solid line 
stands for the ICPA results while the dotted line stands 
for the CPA. In Fig.6(b), the sohd fine is the ICPA re- 
sults and the dotted line is the experimental curve. The 
CPA results suggest a split-band behaviour in the middle 
of the band clearly separating the Pt contribution in the 
low frequency region from the Ni contribution in the high 
frequency region. The overall contribution from the low 
frequency region is much less than that of the high fre- 
quency region in this system because the low frequency 
region is dominated by the heavier atom Pt which has a 
much lower incoherent scattering length B than Ni, 0.1 
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FIG. 6: (a)Incoherent neutron scattering structure factor 
vs. frequency calculated in the ICPA(solid line) and in the 
CPA(dotted line) in NisoPtso (b) Same plot as (a) in the 
ICPA(soIid line) and experimental results U] (dotted line). 



in comparison to 4.5 for Ni. Including only mass fluctu- 
ations and ignoring the Ni-Pt correlated motion induced 
by the environmental disorder gives rise to this spurious 
gap in the CPA results. This point is further discussed 
below in connection with our analysis of the coherent 
scattering results. On the other hand, by incorporating 
the force-constant disorder, as is done in the ICPA, we 
get rid of this spurious gap and obtain good agreement 
with the experimental results, including the position of 
the right band-edge. The influence of the force-constant 
disorder is demonstrated more prominently in the dis- 
persion curves and the line shapes. In Fig. 7, we compare 
the dispersion curves and widths obtained in the ICPA 
from the coherenet scattering structure factors, using the 
force constants as given above, with those in the CPA, 
using the averages of the same force-constants, and with 
the experimental results [||. The top three panels are 
the results obtained in the ICPA for the three symmetry 
directions. The procedure has already been discussed in 
the previous section. The bottom panels are the CPA re- 
sults. The circles and error bars in the left panels are the 
experimental frequencies and the widths (FWHMs). The 
filled circles are experimental results with more accuracy 
and greater resolution. The shaded regions in the top 
panels and in the leftmost panel in the bottom span the 
calculated FWHM's. The FWHM's for the middle and 
the right panels in the bottom are indicated by the thin 
dotted fines. For the [C,0, 0] and the [C, d C] directions, 
the solid lines represent the longitudinal modes and the 
dotted lines the transverse modes. For the C, 0] direc- 
tion, the solid lines show the longitudinal mode while the 
long-dashed and the dot-dashed curves stand for the Ti 
and the T2 transverse modes, respectively. The experi- 
mental results are available only for the [C, 0, 0] directions 
in this system. The ICPA agrees much better with the 
experiments than the CPA for both the longitudinal and 
the transverse branches. The CPA frequencies are gener- 
ally below the experimental ones at low frequencies and 
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FIG. 7: The solid lines are the L branch in all the three 
panels, the dashed lines are the T branch in the left and the 
right panels. In the central column, the long-dashed curves 
are the Ti branches while the dot-dashed curves are the T2 
branches. The shaded regions span the FWHMs. The cir- 
cles in the left panels are the experimental data The 
filled ones are those with better resolution and accuracy. (Top 
panel) Dispersion curves for NisoPtso calculated in the ICPA. 
(Bottom panel)Dispersion curves for NisoPtso calculated in 
the CPA. Here, the shaded regions in the left panel span the 
FWHMs. In other two panels the thin dotted lines denote the 
FWHMs. 



above the experimental ones at high frequencies. The 
discrepancy gets worse as we move from the middle of 
the zone towards the zone edge. This can be explained 
the following way: the high wave-vector region is domi- 
nated by the lighter atoms. The use of the average force- 
constants by the CPA coupled with the domination of 
the lighter mass pushes the frequencies further up, thus 
producing a significant deviation from the experimental 
ones. The severity of this effect can be understood from 
the widths as well. In the CPA, the experimental points 
stay well outside the disorder-induced widths centered at 
the peak frequencies. The discrepancy is substantially 
reduced by the inclusion of force-constant disorder, as is 
seen from the ICPA results. Its inclusion changes the 
dispersion curves qualitatively as well. In the CPA, the 
bands extend fully across the Brilluoin zone for all sym- 
metry directions while in the ICPA, the Pt-dominated 
peaks disappear at high-^ and the Ni-dominated peaks 
wash out at low-(^ for all modes and symmetry direc- 
tions, an effect observed in the experiments. This is a 
clear consequence of the force-constant disorder which 
can be understood by inspecting the spectral line shapes. 
In Fig. 8 (a) we present the partial and the total spectral 
densities for three different wave- vectors, one on the lower 
side of the zone, one in the middle and one at the bound- 
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ary. The solid lines are the total spectral functions, the 
dotted lines are the Ni-Ni partial spectral functions, the 
long-dashed lines are the Pt-Pt contributions while the 
dot-dashed lines are the Ni-Pt contributions. In Fig. 8(b) 
we present the partial and the total coherent scatter- 
ing structure factors i.e. the spectral functions weighted 
by the coherent scattering lengths of the species. The 
coherent scattering lengths of Ni and Pt differ by only 
7%(the scattering length of Ni is 1.03 while that of Pt is 
0.95), much less than do those of Ni and Pd. However, 
even this small difference produces significant changes in 
the line shapes and in the peak frequencies. In Fig. 8(b), 
the solid lines give the total structure factor, the dot- 
ted lines are for weighted Ni-Ni contributions, the long- 
dashed lines give the weighted Pt-Pt contributions and 
the dot-dashed lines stand for weighted Ni-Pt contribu- 
tions. A close inspection of the various contributions re- 
veals the fact that unlike in NiPd, the Ni-Pt contribution 
plays the key role in determining the weight in the middle 
of the band (and in obtaining the merged bands in Fig. 6) 
as well as adding or subtracting weights to the Ni-Ni or 
Pt-Pt contributions, thus elevating or suppressing one 
of the peaks. For example, in Fig. 8(a), in the [.5,0,0]T 
curves, the Ni-Pt contribution adds weight to the total 
spectral function on top of the Pt-Pt peak at the low 
frequenies while it subtracts weight from the Ni-Ni con- 
tribution at higher frequencies thereby causing a weakly 
defined peak at high frequencies. In the [.5,0,0]L and in 
the [1,0, 0]T curves, the Ni-Pt contribution adds weight 
between the Pt-Pt and Ni-Ni peaks, thereby removing the 
gap in the spectrum. The Ni-Pt contribution is totally 
due to inclusion of force-constant disorder, since diago- 
nal disorder produces no such contribution. Thus, the 
CPA produces spectral functions having two well-defined 
peaks corresponding to the Ni-Ni and the Pt-Pt contribu- 
tions with a gap in between resulting in extended disper- 
sion curves and split bands. The effect of incorporation 
of the difference in scattering lengths can be seen from 
these two figures as well. For example, in the [1,0, 0]T 
curves, there are two well-defined peaks in the total spec- 
tral functions, whereas the low- frequency peak is trans- 
formed into a shoulder in the total structure factor. This 
is because Ni has the larger scattering length which en- 
hances the weight associated with the Ni-Pt contribution 
thereby cancelling more effectively the contribution from 
the Pt-Pt part. Similar effects are seen in the [.5,0,0]L 
and [1,0, 0]L curves. Moreover, this weighting sometimes 
produces a weakly defined peak whose FWHM cannot 
be well determined, which explains the observed washing 
out of the dispersion curves noted above. The effect of 
the small difference in scattering lengths is amplified by 
the force-constant disorder through the Ni-Pt structure 
factor. 

In sum, the force-constant disorder plays a significant 
role in Ni5oPt5o, and a theory with mass disorder only, 
like the CPA, fails both qualitatively and quantitatively 
in such cases. On the other hand, the ICPA successfully 
explains the effects of force-constant disorder through its 
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FIG. 8: (a)Partial and total spectral functions calculated 
in the ICPA for various C values in the [C,0,0] directions in 
NisoPtso. (b) Partial and total structure factors calculated 
in the ICPA for various values in the [C,0,0] directions in 
NisoPtso. The solid lines are the total contributions, the dot- 
ted lines are the Ni-Ni contributions, the long-dashed lines 
are the Pt-Pt contributions and the dot-dashed lines are the 
Ni-Pt contributions. The details are given in the text. 



We have presented a straight-forward and tractable 
formulation of the KLGD |l^ method for single-site scat- 
tering of phonons in three dimensional lattices. We have 
demonstrated how this multiple-scattering based formal- 
ism captures the effects of off-diagonal and environmental 
disorder. The use of augmented-space to keep track of 
the configurations of the system has made the formalism 
simple yet powerful. The resulting translational invari- 
ance makes it numerically tractable as well. A signifi- 
cant contribution beyond Ref.18 is the derivation of the 
partial Green's functions in real-space and the deriva- 
tion of the partial spectral functions as well as their sum 
rules. This enables one to make direct comparison with 
neutron scattering experiments because of the incorpo- 
ration of the scattering lengths of the different species. 
We have applied the formalism to real random alloys for 
the first time. In Ni55Pd45 we have demonstrated that 
mass disorder plays the prominent role, and the CPA 
consequently does a rather good job whereas the mean- 
crystal model requires erroneous fitted force-constants. 
Our partial structure factors enable us to understand the 
insensitivity of the normal modes towards the difference 
in the coherent scattering lengths of the two species de- 
spite the significant difference of 43% in this system. The 
NisoPtso results demonstrate the prominence of force- 
constant disorder even in a case where the mass ratio 
is 3. The ICPA agrees well with both the coherent and 
the incoherent scattering experiments, whereas the CPA 
fails, both quantitatively and qualitatively. We are able 
to establish the role of the force-constant differences be- 
tween the species in great detail with the help of the 
partial spectral functions and the partial structure fac- 
tors. We have consequently clearly demonstrated that for 
systems like NiPt, where the force-constants are strongly 
species dependent, the determination of their values is 
crucial. However, we had no prior information about 
the species dependence of the force-constants. Intuitive 
arguments led to a set of force-constants which turned 
out to be quite good. A better understanding of the 
role of disorder in the lattice dynamics of random al- 
loys could be achieved with prior information about the 
force-constants. These could be obtained, e.g., from first 
principles calculations on a set of ordered alloys. 
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effect on the partial structure factors. It also demon- 
strates the relative importance of the contributions of 
various atomic species to the coherent and incoherent 
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APPENDIX A: MATRIX ELEMENTS 

For the calculations in the nearest-neighbor approxi- 
mation, one needs to evaluate the {3 {Z + 1)}^ matrix 
elements of the operators K', K'^, V (g) and G"^^^ and 
use them as inputs. These evaluations are done in aug- 
mented space using Eq.(14). The symmetry of the lattice 
structure is used to reduce the number of matrix elements 
evaluated. Here, we give results only for an fee lattice. 
All the matrices are, therefore, of dimension 39 x 39. 

In an fee system, each atom has 12 nearest neigh- 
bors with coordinates (±|,±5,0), (±^,0,±5), and 
(O, ±i,±i) with respect to the coordinates of the ref- 
erence atom at (0,0,0). The force-constants, between 
the atom and its neighbors satisiy the following cubic 
symmetry relation. 

^o/ = <" = < = 4t> (Al) 



where i?oj = -^Rok and k and j are two neighbors on 
opposite sides of site 0. For example, the force-constant 
matrix between the atoms (0, 0, 0) and (5,5,0) is of the 
form, 



/a b 0\ 

(A(ooo iio) = { b a \ . (A2) 
\0 g J 



The force-constant matrices between the atom and its 
other neighbors can easily be calculated from (A2) via 
the cubic symmetry operations. The results are: 



a/3 
ij 



g^xx 



rrvjj 
= 
= Df ^ 
= Df" 

= 2R^ X 2R^ X 
= 



4^/xx 



iU=j, a = (3, 
iii = j, /3, 

if i = 0, j = n = 1 - 12, a = (3, Rf = 0, 
if i = 0, j = n = 1 - 12, a = P, Rf , 

ifi = 0,j = n=l-12,a^P, 
otherwise, 



(A3) 



and 



= {GvD'^i = (^^ca)f/ ' ah i,j,a, and (3. 



Also, we find 



= 



^2 

2Rf X 2i?f X 

3 3 ^ 



= 



a/3 



if i = j = 0, a = /3, 
if i = j = 0, a ^ /3, 

if i = 0, j = n = 1 - 12, a = /3, ii? = 0, 
if i = 0, i = n = 1 - 12, a = /3, ii^ ^ , 
if i = 0, j = n = 1 - 12, a ^ /J, 
if i ^ 0, j = i = n = 1 - 12, 
otherwise. 



(A4) 



and 



a/3 



a/3 



for all i,j,a, and /?. 
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Similarly, we obtain 

= 4D^ySin{q^)Sin{q0) if i = j = 0, a ^ 

= D5 if i = 0, j = n = 1 - 12, a = /3, i?^ = 0, 

= De if i = 0, j = n = 1 - 12, a = /?, i?^ ^ , 

= 2Rf x2R'^ xDr if i = 0, j = n = 1 - 12, a ^ f3, 

= Df ^ - £>f ^ if i ^ 0, i = i = n = 1 - 12, a = /?, i?," = 0, 

= _ ifijLO, j = i = n=l-12,a = /3, Rf j^O , 

= 2R^ X 2iif X {Dl^ - D=[y) iii^O, j = i = n= 1-12, (3, 

= -Df^e'^-^' a j ^ i = n = 1 - 12, Rj = -Ri, a = f3, Rf = 0, 

= _£)|^e'«-^' if j ^ i = n = 1 - 12, = -^i, a = p, i?^ ^ 0, 

= -Dfe'^-^' X 2i?f X 2i?f if j ^ i = n = 1 - 12, Rj = -Ri,aj^ (3, 

= otherwise, (A5) 



and 



V {Q)tf = (9)*) = y ®? , for all i,j,a, and /3. 



In these evaluations, we have used the notation Rqj = Rj, n = 1 — 12 to represent the 12 nearest neighbors, and the 
notations 



m 




A , B 

CAfn + cbiti , 




m' 




y/cACB {m^ ~ m^) , 




rh 




cbiti^ + CAin^ , 




J^XX 




c\aAA + c\aBB + 2cACBaAB, 








CaQAA + C%gBB + 2CACBgAB, 








c\hAA + C%hBB + 2CACBbAB-, 








\/CACB{cAaAA - CBaBB + {cb 


- CA)aAB}, 


J^IXX 




y/CACB{cAgAA - CB9BB + {cb 


- CA)gAB}, 






y/CACB{cAbAA - CBbBB + {cb 


- CA)bAB}, 


T~iXX 




caCb {aAA + aBB - 2a ab) , 




^3 




caCb {9aa + Qbb - 2gAB) , 








CaCb {bAA + bBB - 2bAB) , 




T^XX 




CACB (ttAA + CLbb) + {c\ + C%) 


aAB, 






CACB [gAA + gBB) + {c\ + c|) 


9AB, 






CACB {bAA + bBB) + {ca + c|) 


D5 












- ) + i?^^e*«•^^ 




D7 












(m - m) - 8 (D^^ - Df") - 


4 (Df ^ - d: 



-ADf^'Cosq^Cosqs; 1,^1^01. (A6) 



The symmetries of the force-constant matrices are re- is captured in F (g)^'^ through the quantities Dg . When 
fleeted in the operators as well. The effect of itineration there is no force-constant disorder the terms vanish. 
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and V becomes independent of q. A (f-indcpendcnt self- Using Eqs.(6)and (11), 
energy results, and we arrive at the CPA equations. To 
illustrate how to obtain the various matrix elements of 
the operators, we present the calculation of where 

iii = (i,i,0). (0/|m|l/o) = 0, 

= (0/|K|l/o), 

= (0/1 {muP' - |l/o). and using Eqs.(7) and (11), 



(0/|$|l/o) = (/|$oi|/o), 

= ((V^(Ao| + ^{B^\) , (V^(^i| + V^(i?i|) \{<\>^^n^n^ + 0oT^?(f 5?f + '/'^i^^^^yf + 

(V^l^o) - V^l^o)) , (V^l^i) + V^l^i))). 



If we use the cartesian coordinates explicitly, then the xx 
component is, for example, 

(/l*Oll/o) = \fcAC^{CAaAA- CBaBB^{CB- Cjs)aAB) 

= Df". 

The other components can be calculated similarly. 
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